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We show that a de Sitter spacetime is a solution of Einstein's field equations with the energy 
momentum tensor of a self-interacting, classical Maxwell-Boltzmann gas in collisional equilibrium. 
The self-interaction is described by a four-force which is quadratic in the (spatially projected) particle 
four-momenta. This force does not preserve the particle number and gives rise to an exponential 
increase in the comoving entropy of the universe while the temperature of the latter remains constant. 
These properties of a gas universe are related to the existence of a "projector-conformal" timelike 
Killing vector representing a symmetry which is "in between" the symmetries characterized by a 
Killing vector and those characterized by a conformal Killing vector. 
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I. INTRODUCTION 

Standard inflationary cosmology relies on the dynam- 
ics of a scalar field with a suitably designed potential term 
allowing for an effectively negative pressure of the cos- 
mic substratum Alternatively, there have been nu- 
merous attempts to describe early phases of accelerated 
expansion, either power-law or exponential, as nonequi- 
librium (imperfect) fluid phenomena (see |||| and ref- 
erences therein) . Studies along this line make use of the 
fact that fluid viscosities are dynamically equivalent to ef- 
fective negative pressures. Such kind of pressures occur, 
e.g., due to internal interactions inside relativistic gases 
However, the dissipative fluid approaches face the 
following general problem. Their reliability is restricted 
to small deviations from perfect fluid behaviour, i.e., 
close to (local) equilibrium, while inflation in this context 
is necessarily a far-from-equilibrium phenomenon. Al- 
though the imperfect fluid dynamics was shown to admit 
inflationary solutions (see |],|| and references therein) 
their physical significance is unclear ||. The question 
whether or not bulk viscosity may drive inflation || re- 
mains open. 

A different line of fluid dynamical early universe inves- 
tigations implying effective negative pressures as well is 
connected with cosmological particle production jio] . 
Here, the negative bulk pressure is not due to devia- 
tions from (collisional) equilibrium but is a consequence 
of the phase space enlargement of the fluid particle sys- 
tem. Negative pressure terms of this kind in the lo- 
cal energy-momentum conservation relations may be re- 
garded as equivalent to corresponding not necessarily 
small source terms in the energy-momentum balance of a 



perfect fluid. In other words, the effective nonequilibrium 
description in terms of imperfect fluid quantities mimics 
a perfect fluid with continuously increasing particle num- 
ber, supposedly of quantum origin. While the magnitude 
of conventional viscous pressures which are due to par- 
ticle number preserving interactions within the fluid is 
severely restricted by close-to-equilibrium conditions, the 
magnitude of effective, cosmological particle production 
describing bulk pressures is not. 

Previous investigations have shown that such kind of 
(generalized) equilibrium particle production (creation of 
perfect fluid particles with minimal entropy production) 
may substantially modify the standard cosmological dy- 
namics, including the possibility of "reheating" phenom- 
ena and power-law inflation ^(j , g2| -g4[. Particularly in- 
teresting connections were established between spacetime 
symmetries, described by a conformal timelike Killing- 
vector and the production rate for particles at equilib- 
rium. Depending on the equations of state the symmetry 
requirements turned out to fix the creation rate for fluid 
particles at equilibrium. Moreover, for a gas universe the 
production of particles could be traced back to specific 
self-interacting forces on the (classical) microscopic con- 
stituents of the cosmic medium p4fl . Essential features 
of the cosmological dynamics may be discussed in terms 
of microscopic particle motion in a (classical) force field. 

In section III of the present paper we identify those 
forces on the particles of a simple relativistic gas gov- 
erned by an equilibrium distribution function which, on 
the phenomenological level, give rise to an effective bulk 
pressure 7r = — (p + p), where p is the energy density of 
the cosmic medium and p is the corresponding equilib- 
rium pressure. A negative pressure 7r of this magnitude 
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is equivalent to a particle production rate that coincides 
with the expansion rate of the universe (see Eq. be- 
low). We show that a gas configuration such as this is 
intimately connected with a symmetry of the spacetime 
characterized by a "projector-conformal timclike Killing- 
vector" (PCTKV) £ a which we generally define by 



£ c „9ik = Ci;k + Ck;i = 2<j>P$ 



(1) 



where £ ( gi k is the Lie-derivative of the metric tensor gi k 
with respect to Co and <j) = 4> (x) is a spacetime function. 
The quantity Pi k is the projector 



Pik = 



lik 



CaC a 



P lk ( k = o 



(2) 



on surfaces orthogonal to the timelike vector ( a . The spe- 
cial case of interest here corresponds to Ca = u a /T, where 
u a is the macroscopic four- velocity of the cosmic medium 
and T is its equilibrium temperature. The projector Pi k 
then reduces to hi k = gtk + UiU k , the projection tensor 
orthogonal to u l . The first part of the paper (section 
II) therefore intends to clarify the general implications 
of the PCTKV property of u a /T for the relativistic fluid 
dynamics. Section IV considers the case that a medium 
such as characterized in section III dominates the dynam- 
ics of a homogeneous and isotropic universe, while section 
V summarizes our main findings. Units have been chosen 
so that c — ks = Ti = 1. 



II. FLUID DYNAMICS WITH A 
"PROJECTOR-CONFORMAL" TIMELIKE 
KILLING- VECTOR (PCTKV) 

In fluid spacetimes the Lie derivative £^g a b of the 

metric tensor g ab with respect to u a /T, where u a u a = 
— 1, plays a well-known role to characterize symmetries 
of the metric. Especially interesting are cases in which 
u a /T is cither a Killing vector or a conformal Killing 
vector. Like any symmetric tensor, the Lie derivative 
of the metric may generally be split into contributions 
parallel and perpendicular to the four- velocity: 

with 

C ik = 2Au.Uk + B t u k + B k Ui + 24>h ik + b lk (4) 

and hi k u k = Biu 1 — bi k u k = b\ = 0. Combining relations 
(§) and (0) we have 



A = \f lk u l u k = , 



-h l m u Ca 



1 
T 



T 



(5) 



(G) 



and 



b a b = 



h\h\ - -h ab K k 



a 



ik 



1 

T 



V b u a + W a u b - -@h a b 



(8) 



where & = u\ is the fluid expansion, T = T a u a , V a T = 
/i a Tft, ii a is the four-acceleration and a a b is the shear. 

According to different choices for the quantities A, B m , 
(j), and b ab we introduce the following classification for 
Ua/T: 

(i) u a /T is a Killing- vector for A = B m = = b ab = 0, 
corresponding to 



TJ:k VT ' 



(9) 



6 = T = 



T 



u a = , a ab = . (10) 



J tn u ab 



(ii) u a /T is a conformal Killing-vector for A = —(f) and 

(ii) 



equivalent to relations (Q) and 

f e v a r 



T 



T 



Ua = , (Tab = 



(12) 



(iii) u a /T is a "projector-conformal" timelike Killing- 
vector (PCTKV) for A — B m = b a b = 0, i.e., 



or relations (^) and 



T 

- = , 
T T 



U a = , cr ah = 



(13) 



(14) 



In this paper we are interested in the third of these 
cases. In the following section we will give a microscopic 
realization of the PCTKV behavior on the level of kinetic 
theory and identify a spacetime which admits a PCTKV. 

In the present section we show that the property (j?]) 
with a nonvanishing together with Eqs. ( |l4| ) may be 
realized in case the fluid particle number is not preserved. 

Introducing a length scale a according to = 3a/ a, 
the number N of particles in a comoving volume a 3 is 
N = rta 3 , where n is the particle number density. The 
corresponding particle number flow vector is N a = nu a . 
Denoting the phcnomenological particle production rate 
by r, the particle number balance may be written as 



2 



iV 

N a =n + On = n—= nT 
N 



(15) 



The energy-momentum tensor of a general fluid is 
given by 



T ik _ T ik + nh ik + n ik + q i u k + q k ul 



with T& 



pu x v^ + ph lk and w^Uk 



qui 



h lk Ui = 0. Local energy-momentum conservation T. k = 
implies 

p = -6 (p + p + tt) - V Q <7 a - 2u a9 a - a Qb ^ ab . (17) 
With the help of the Gibbs equation (see, e.g., |2jJ) 

(18) 

where s is the entropy per particle, we obtain 



Tds = d— + pd— . 

n n 



nTs = p — (p + p) 



(19) 



Taking into account the standard definition S l — nsu 



q /T of the entropy flow vector S % and using Eqs. ( |15[ ) 
@, and @ yields 



P+P, 



T 



rjilh rjiih 

1 ~ J (o) 



Ui 

TJ;k 



(20) 



Under the PCTKV conditions (13h, Eq. (EOh reduces to 



nsr 



P + P 
T 



r 



(21) 



For q a — tt = the right-hand side of this equa- 
tion coincides with ns. The especially interesting case 
07r = — (p + p) r which makes the right-hand side of 
Eq. (|Tj) vanish implies = nsT (and s = for 

q a = TT ab = 0): There is entropy production only due to 
the enlargement of the phase space of the system but not 
to conventional dissipative processes within the fluid, i.e., 
the entropy production is minimal. States in which the 
right-hand side of Eq. ( pil| ) vanishes (and T > holds) 
are states of "generalized equilibrium" |23j . This kind of 
equilibrium was originally introduced on the basis of the 
conformal Killing vector property ([ll]) of Ui/T. Here we 
enlarge the generalized equilibrium concept to include 
the PCTKV case ( |l3| ) as well. We consider a fluid to 
be in "generalized equilibrium" if (i) S l i — nsT = with 
r > holds and (ii) u%jT satisfies either the conditions 
(li), equivalent to Eqs. nA) and (UW, or the conditions 
(Jq), equivalent to Eqs. (fx) and (llj). 

For massive particles in a homogeneous and isotropic 
universe the first case, minimal entro py p roduction un- 
der the conditions (RW , dealt with in fe3U24f| , was shown 



to imply power-law inflation. In this paper we are in- 
terested in the second case, i.e., minimal entropy pro- 
duction according to S\ — nsT = together with the 
PCTKV property (|l|).' We will demonstrate that this 
type of generalized equilibrium under the conditions of 
homogeneity and isotropy requires a de Sitter universe 
for arbitrary fluid equations of state. 

It is a characteristic feature of generalized equilibrium 
of both types that the particle production rate is not 
an arbitrary parameter but determined by consistency 
requirements. Given equations of state in the general 
form 



p = p T) 



p = p(n,T) , 



(22) 



differentiation of the latter relation and using the bal- 
ances (^5j) and (|l7]) yields 



T 
T 



dp/dT ' 



(23) 



where the abbreviations 
(d P /dT) n 



dp 

dp 



(dp/dT) n 



dp 
dT 



dp_ 

dT 



have been used. Restricting ourselves to the case s = 0, 
the temperature laws in Eqs. ( p^ ) and (^3|) are only 
consistent for T = 6. Via the relation 



9tt = -(p + p)T 

which makes the right-hand side of Eq. 
viscous pressure tt is fixed as well: 



r = e 



TT = - (p+p) 



(24) 
vanish, the 

(25) 



From Eq. (|15|) follows that n vanishes. Because of the 



equations of state (22), constant values of n and T along 



the fluid flow lines imply constant pressure and constant 
energy density as well: 



- = £ = ° 
n T 



=> p = p = 



(26) 



We recall [g3[ that generalized equilibrium under the 
conformal Killing- vector conditions (11) is characterized 
by Eq. (pi with 



instead of Eq. fl25|). 

Statements on the spatial dependences of the thermo- 
dynamic quantities may be obtained from the momentum 
balance 

(p + p + n) it m + V m (p + tt) = , (27) 

where we restricted ourselves again to q a — TT ab = 0. 
Because of relation (pa) we find 
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-V„7T 



V a p = 



(28) 



The energy density is also spatially constant, the pressure 
not necessarily. With the help of the second equation of 
state (22) the condition ( p8| ) provides us with a relation 
between V a T and V a n: 



V a T n dp/dn V a rt 
T ~ ~Td~pJdf~n~ 



(29) 



where 



dp 

dn 



p + p 



T dp 
ndf ' 



which is a consequence of the fact that the entropy is a 
state function. 

Restricting ourselves to a classical gas with p = nT 
and using the first relation ( p8| ) and Eq. ( f29f ) as well as 
the second relation ( |l4|) we obtain 



nT 



Tdp/dT _ 
n dp/dn 
T dp/dT 



V T 
T 



n dp/dn 

From the Gibbs-Duhem relation 
dT 



dp=(p + p) — +nTd(^) 



follows under such circumstances that 



Tdp/dT 
n dp/dn 



V T 
T 



(30) 



(31) 



(32) 



The Gibbs equation (|T|) together with the second re- 
lation (J2q) then provides us with 



V„* = - ~ + l 



= - ^ + 1 



n 

T dp/dT 
n ndp/dn 



(33) 



The spatial gradient of the entropy per particle does not 
vanish in general. 

We finish this section by considering the general rela- 
tions (pq ) - (|33|) for the limiting cases of pure radiation 
(ultrarelativistic matter) and nonrelativistic matter. 

(i) p — nT, p = 3nT (ultrarelativistic matter): 
The relation V m p = implies V m T/T = —\7 m n/n, 
equivalent to V m p — — V m 7r = 0. From Eq. d32| ) one 
obtains 



P= 3 



(34) 



Although the spatial pressure gradient vanishes, the cor- 
responding gradients of the temperature, the number 
density and the entropy per particle are different from 
zero unless the fluid motion is geodesic. 



(ii) p = nT, p = nm 
matter) : 

In this case the second relation 



T 

and we find 



2 m V m n 

3 T n 



nT, m 3> T (nonrelativistic 
reduces to 

(35) 



(m»T) , 



V TO p = -V m 7T 



-p u„ 



(m»T) 



(36) 



Nonrelativistic matter (p 3> p) in the present case com- 
bines a homogeneous energy density with a generally in- 
homogeneous equilibrium pressure. Only the total effec- 
tive pressure is homogeneous as well. 
For the spatial gradient of p/T we get 



m V m T _ 
T T 
' p — m 
T 



m 



, (m > T) . (37) 



It is well known that p — m is the nonrelativistic chemical 
potential (see, e.g., p6[ ). 

In the following section we show how a gas character- 
ized by the properties ( p4| ) and fl25| ) may be realized with 
the help of relativistic kinetic theory. 



III. KINETIC THEORY FOR A GAS IN A FORCE 
FIELD 

A. General relations 

The conventional kinetic theory of a simple relativistic 
gas relies on the concept of pointlike particles which may 
interact through elastic, binary collisions. Inbetween the 
collisions which are assumed to establish an (approxi- 
mate) local or global equilibrium of the system the par- 
ticles move on geodesies of either a given spacetime or 
a spacetime selfconsistently determined by the gas parti- 
cles themselves. Sophisticated solution techniques for the 
corresponding Boltzmann equation have been developed 
and applied to numerous physically relevant situations 
p5| p9[ . The usual procedure here is first to character- 
ize equilibrium states, i.e., states with vanishing entropy 
production and to relate the parameters of the corre- 
sponding distribution function to macroscopic (perfect) 
fluid quantities. Noncquilibrium situations are then, in a 
second step, taken into account as deviations from equi- 
librium. Obviously, geodesic particle motion is a highly 
idealized case. In reality, particle worldlines are supposed 
to deviate from geodesies since the particles will be sub- 
ject to additional interactions in general. We assume 
here that these interactions may be modelled as effective 
forces on the particles. The kinetic theory for particles 
under the influence of various forces was considered, e.g., 
in HH^Q. Following the lines of |H) we will fo- 
cus here on equilibrium states of the gas in a force field. 
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Our main objective will be the characterization of equi- 
librium configurations of a gas under the influence of a 
self-interacting force which is quadratic in the particle 
four-momenta. By "self-interacting" we mean that the 
force, except its dependence on the microscopic particle 
momenta, also depends on macroscopic quantities, char- 
acterizing the gas system as a whole (see below). This 
force will neither preserve the particle number nor the 
energy momentum and it will give rise to entropy pro- 
duction. In particular, it will turn out that a quadrati- 
cally self-interacting force realizes states of "generalized 
equilibrium" characterized below Eq. ( pl| ) , including the 
condition s = for "adiabatic" (or "isentropic" ) particle 
production. 

The one-particle distribution function f — f (%,p) for 
relativistic gas particles under the influence of a four- 
force F 1 = F l (x,p) obeys the Boltzmann equation 



. +mF 1 ^- = C[f] , (38) 
<V dp 1 11 y ' 



where / (x, p) p k nkdHdP is the number of particles whose 
world lines intersect the hypersurface element rikdY, 
around x, having four- momenta in the range (p,p + dp). 
The quantity dP — A(j>)6 (p l Pi + to 2 ) dP± is the volume 
element on the mass shell p l pi = —to 2 in the momentum 
space. A(p) — 2, if p l is future directed and A{p) = 
otherwise; dP^ = ^/^gdp a dp 1 dp 2 dp 3 . 
C[f] is Boltzmann's collision term. Its specific structure 
discussed e.g. by Ehlers |^6| will not be relevant for our 
considerations. Following Israel and Stewart |^7j we shall 
only require that (i) C be a local function of the distribu- 
tion function, i.e., independent of derivatives of /, (ii) C 
be consistent with conservation of four-momentum and 
number of particles, and (iii) C yield a nonnegative ex- 
pression for the entropy production and do not vanish 
unless / has the form of a local equilibrium distribution 
function (see ( [iij ) below). Equation (38) implies that the 
mass-m particles inbetween the collisions move according 
to the equations of motion 



dx % 
d7 



P 



Bp 1 
dj 



= mF l 



(39) 



where 7 is a parameter along their worldline which for 
massive particles may be related to the proper time r by 
7 = t/to. Since the particle four-momenta are normal- 
ized according to p l pi = —to 2 , the force F l has to satisfy 
the relation piF 1 = 0. 

Both the collision integral C and the force F % describe 
interactions within the many-particle system. While C 
conventionally accounts for elastic binary collisions, we 
intend F % to model different kinds of interactions in a 
simple manner. Strictly speaking, F 1 should be calcu- 
lated from the microscopic particle dynamics and, conse- 
quently, depend on the entire set of particle coordinates 
and momenta characterizing the system of gas particles. 
We will introduce here the simplifying assumption that 



F z be an effective one-particle quantity which instead of 
depending on the coordinates and momenta of the re- 
maining particles, is supposed to depend on macroscopic 
fluid quantities characterizing the system as a whole. At 
the moment we do not specify this force. It will be de- 
termined below by general equilibrium conditions. We 
expect the concept of a self-interacting force to be useful 
in circumventing some of the general problems inherent in 
attempts to formulate a relativistic statistics for interact- 
ing many-particle systems (see, e.g., ||^] and references 
therein) . 

The particle number flow 4-vector N l and the energy 
momentum tensor T lk are defined in a standard way (see, 
e.g., {§) as 



N* = J dP P *f(x,p) , f ik = J dPp l p k f{x,p) 



(40) 



While it will turn out that in the equilibrium case the 
first moment of / in Eq. ( [io| ) may be identified with 
the quantity N l — nu % introduced below Eq. (|l4|), we 
have used here the same symbol immediately. The sec- 
ond moment of / denoted by T lk in Eq. (40) will not, 
however, coincide with the energy-momentum tensor T lk 
in Eq. (|l6|). The integrals in the definitions ( [4b"| ) and 
in the following are integrals over the entire mass shell 
p l Pi = —m 2 . The entropy flow vector S a is given by p6[, 



S a = 



J P a [/ln/-/]dP 



(41) 



where we have restricted ourselves to the case of classical 
Maxwell-Boltzmann particles. 

Using well-known general relations (see, e.g., pq] ) we 
find 



iV a = 

x y ;a 



and 



T ak = J P a (cw-mir^jdP, 

S -a = -j^f{c[f}-mF i ^C) dP 



(42) 



(43) 



In the follo wing we will focus on the force terms in the 
expressions ( |42| ) and (|43|). In order to separate the cor- 
responding contributions from those due to the collision 
integral we will restrict ourselves to collisional equilib- 
rium from now on. Under this condition In / in Eq. ( |4"3| ) 
is a linear combination of the collision invariants 1 and p a 
and the contributions due to C [f] in formulas (42) and 
( [43| ) vanish. The corresponding equilibrium distribution 
function becomes (see, e.g., Pq| ) 



f (x, P )=exp [a + (3 aP a ] 



(44) 
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where a = a (x) and p a (x) is timelike. 

Inserting the equilibrium distribution function ( fli]) into 

Eq.(|38|) one obtains 

P^.a + f3(a;b)P a p b = . (45) 



For a vanishing force F l the latter condition reduces to 
the "global" equilibrium condition of standard relativis- 
tic kinetic theory. For F % ^ condition (45) is a "gener- 
alized" equilibrium condition (see [23 and below). 
Use of Eq. (J44[) in the balances (42) yields 



K a = -mft J F'fdP , 

T"$ = -Tufa J p a F l fdP . 

For the entropy production density (^) we find 

mfc J[a + f3 aP a ]F l InfdP 
~aN a a - [3 a f ab . 



(46) 



With / replaced by f° in the definitions (|40|) and (41), 
N a , T ab and S a may be split with respect to the unique 
4-velocity u a according to 



iV a = nu a , T ab = pu a u b + ph ab , S a = 



(48) 



where u a , h ab , n, p, p and s may be identified with the 
corresponding quantities of the previous section. The 
exact integral expressions for n, p, p and s are given by 
the formulae (177) - (180) in @. 

Using the structure (|48|) for N a and defining 



A / F7°dP , 



(49) 



the first Eq.(|46|) becomes n + Qn = nT [cf. Eq. (|15[)] 
which justifies definition (^|). Similarly with the de- 
composition (E8J) and the abbreviation 



t a = m/3i / p a F l fdP 



we obtain 



implying 



f ab b + t a 



p + e{p + p) =u a t a , 

(p +p)ll a + V a P = -h ai t % 



(50) 



(51) 



(52) 



The energy-momentum tensor T lk is not conserved. Ob- 
viously, the balances (52) following from Eq. (|5l]) are 
identical to the balances 



p + Q(p + p + TT) = 0, 
{p + p + ir)u a + V a {p + ir) =0 



(53) 



following from the local conservation T ab h = of an effec- 
tive energy- momentum tensor T ab , 

T ab = pu a u b + + tt) h ab (54) 

with the identifications 



U a t a 



-Gtt 



hnit 



(55) 



This mapping of the "source" term t a onto an effec- 
tive viscous pressure tt of a locally conserved energy- 
momentum tensor T ab was explicitly shown to be con- 
sistent for specific "sources" t a , depending on the first 
and second moments of the distribution function only 
|]l9| , p0|j23[ . In the following subsection we will show that 
this interpretation continues to hold if specific third mo- 
ments are involved. 

Let us now decompose the four-momenta p a into p a = 
Eu a + \e a where e a is a unit spatial vector, i.e., e a e a = 1, 
e a u a = 0. Consequently, one has E = —u a p a and 
A = e a p a and the mass shell condition p a p a = —m 2 
is equivalent to A 2 = E 2 — m 2 . Moreover, h a \,p a p b = 
(47) _\2 j g va jjj p or force F m we write analogously 
F m = F (x,p) u m + K (x,p) e m where F = -u m F m and 
K = e m F m . The requirement p m F m — will be auto- 
matically fulfilled for 



E 

V 



F(x,p) 



(56) 



With the familiar identification f3 m = u m /T it is obvi- 
ous that only the part u m F m = —F contributes in the 
"sources" ( j49| ) and (|50|). 

B. Quadratic self-interaction 

Any specific force relies on reasonable assumptions 
about F(x,p). In a previous paper we investigated 
the most general linear dependence of F on the particle 
momentum p a . In this paper we assume F to depend 
quadratically on the spatially projected four-momenta A, 



F (x,p) = - Ul F l = XFi (x) + X 2 F 2 (x) 



(57) 



Fi (x) and F 2 (x) are spacetime functions to be deter- 
mined by the equilibrium conditions of the gas. It will 
turn out that it is just this force under the action of which 
the particle energy E is preserved in a homogeneous and 
isotropic universe [cf. Eq. (|^) below]. 
An equivalent way of writing the force (|5^) with F given 
by Eq. © is 

F l (x,p) = g abP a [u l e b - u b e L ] [F x (x) + XF 2 (x)} . (58) 

The equilibrium condition (J45|) becomes 



P a u,a + P { a,b)P a p b = ^ [ e -P a ^ ( X ) + h *bP a p b F 2 (X)] 



(59) 
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It is satisfied for 



and 



a = , e a V a a = —Fi 



P{a-b) = 4>{x)h ab , 



(60) 



with <j) = fr-Fi- Since, on the other hand, relation (0) 
holds we find 



F2 = 

dm 



(62) 



Identifying f3 a with u a /T, the equilibrium condition 
/or a Maxwell-Boltzmann gas in the force field ( p^ ) 
coincides with the PCTKV condition fiT3j). 

We recall that an ansatz for F linear in the particle 
momenta instead of the structure (|57]) analogously repro- 
duces the conformal Killing- vector condition ( p"T| ) p4j . 

Through F 2 the force depends on the fluid expansion, a 
quantity characterizing the g whole on the macro- 

scopic level. Since both the microscopic particle mo- 
menta and macroscopic fluid quantities of the system of 
gas particles enter the four-force, the latter represents a 
self-interaction of the gas. 

Insertingthe force (|58|), equivalent to the expressions 
( |56| ) and (p7|), into the source term ( [49| ) and using 3p — 
h ab f ab = 3nT we obtain 



T = 3mF 2 (x) 



Together with F% from Eq. (|62j) this implies T = 
Eq. (p5|)] and, consequently, the relation h/n = 
For the source term <M\) one obtains 



(63) 

e [cf. 



777 ~ 777 

*° = - y Fi (x) e t T al - -F 2 (x) h m M aim , (64) 

where M aim = J &Pf°p a p l p m is the third moment of the 
equilibrium distribution function. The projected source 
terms in the balances (B3) are 



Uat a = 3n-F 2 (x) , 

where we have introduced the quantity 

V= -\u a h lm M am , 



and 



h na t a = -e n p—F x (x) . 



(65) 



(66) 



(67) 



The third moment M am enters the energy balance but 
not the momentum balance. The momentum balance in 
Eq. ( |52] ) with the "source" term ( p7| ) coincides with the 
corresponding balance in p3| . It follows that the argu- 
ments in p3[ which prove the consistency of the momen- 
tum balances in Eqs. (E^) and (|53|) together with the 



second relation in Eq. (|55|) apply in the present case as 
well. Realizing that 



u r> a f = 

aP 1 



(68) 



(61) (the derivative has to be taken for a = const) , we find 



1 , dT r ' 



dp 



3 " mn d{$) 



Using here for the equilibrium pressure p 
_ 47rm 4 K 2 (f ) 



t kJ- exp \a] 

(2irf (f) 2 



together with 



d (K 2 (z)\ K 3 (z) 



dz V z 2 



(69) 



(70) 



(71) 



{Ki and Kj, are Bessel functions of the second kind) , one 
obtains 



rj = nTm -A^ri = T (p + p) 



(72) 



where we took into account that the enthalpy per par- 
ticle h is given by h = (p + p) jn = mK 3 (^) / K 2 (^). 
Consequently, the source term u a t a becomes 



i a t a = 3m {p + p) F 2 (x) = {p + p)Q . 



(73) 



Together with the first relation in Eq. ( J55[ ) this is con- 
sistent with 7r = — (p + p) [cf. Eq. (Pq)]. 

With the identifications a = p/T, [3 a = Ug/T and 
s = P^jr — j;, the entropy production density (|47|) may 
be expressed in terms of the sources T and u a t a to yield 



S a „ - nsT 



P + P^ , u at a 



T 



T 







(74) 



Recalling that Eq. (|6l]) implies T = O and n = 0, we 
have derived all the generalized equilibrium conditions 
discussed on a phenomenological level in the previous 
section. It is remarkable that the force (^6|) with the 
quadratic dependence ( p7| ) does not only lead to condi- 
tions (^) but automatically guarantees the vanishing of 
S a a — nsT, i.e. s = 0, resulting in T/T = according to 
the temperature law (p3|). XTie condition s = /or 'Wi- 
abatic" or "isentropic" particle production is a property 
of the force considered here and needs not to be postulated 
separately. 

For any T = 6 > Eq. @ implies S. Q a > 0. T/ie 
self-interacting force provides us with a nonnegative ex- 
pression for the entropy production in an expanding uni- 
verse. 

This completes our microscopic derivation of general- 
ized equilibrium characterized by Eqs. ( pl[ ) and ( p4| ) un- 
der the PCTKV condition (O). 
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IV. THE SELF-INTERACTING GAS UNIVERSE 

Having clarified the implications of generalized equilib- 
rium, especially the consequences of relation (|l3|), both 
phenomenologically and on the level of kinetic theory, we 
now consider Einstein's field equations with the energy- 
momentum tensor (^4|). This corresponds to a situation 
where matter in generalized equilibrium dominates the 
dynamics of the universe. Since generalized equilibrium 
was shown to be realized by a self-interacting gas we call 
such kind of configuration a self-interacting gas universe. 
While it is generally an open question to what extent 
the hot and dense early universe is accessible to a ki- 
netic description, a gas is the only system for which the 
correspondence between microscopic variables and phe- 
nomenological fluid quantities is sufficiently well under- 
stood. This makes gas universes interesting toy models 
and we hope that such kind of approach also gives an 
idea of the relevant physics in our real universe. 

The target of this section is to demonstrate explicitly, 
that the specific self- interaction (|5q ) (equivalent to the 
combination of Eqs. (56) and (|57|)) under the equilib- 
rium conditions (^y) and (|6|) allows us to exactly inte- 
grate both the equations for the cosmic scale factor of a 
homogeneous and isotropic self-interacting gas universe 
and the corresponding microscopic equations of motion 
for the individual gas particles. A self-interacting gas 
universe represents an exactly solvable model both mi- 
croscopically and on the phenomenological fluid level. 

In the spatially homogeneous case the function F\ van- 
ishes [cf. Eq. (pQ)]. The lenght scale a coincides with the 
scale factor of the Robertson- Walker metric and obeys 
the equations 



3— = np , 
a A 

where n is Einstein's gravitational constant, and 



(p + p + ir) 



(75) 



(76) 



Together with the expression ( p5| ) for n the last equation 
yields a/a = H = const, where H is the Hubble pa- 
rameter, implying an exponential behaviour of the scale 
factor, a <x exp [Ht], 

A homogeneous and isotropic simple gas universe with 
arbitrary equation of state and quadratic (in the spa- 
tially projected microscopic particle four-momenta) self- 
interaction inbetween elastic binary collisions requires a 
de Sitter spacetime to be in (generalized) equilibrium. 

Evidently, this also implies that the de Sitter met- 
ric admits a PCTKV. The consistency of this statement 
may be checked from Eq. ( |l3| ) directly by using that 
the temperature T is constant both in space and time, 
together with the well-known decomposition of the co- 
variant derivative of the four- velocity |3^J3J . 



Ui-n — UjU Ti 



(77) 



where u> a b = h^h^u^-d]- The present case is characterized 
by u a = o~ a b = ijJab = and = 3^ = const as well as 
(f> = H = const [cf. Eq. (0)]. 

In other words, the selfinteraction of a classical gas is 
able to realize an effective fluid equation of state P e f / = 
p + tt = —p [cf. Eq. (p5|)] which in a cosmological con- 
text is usually obtained with the help of a scalar field. 
A scalar field represents an "exotic" kind of matter with 
mainly theoretical evidence at the present state of knowl- 
edge. We argue here that it may occasionally be helpful 
to have alternative ways of considering issues of inflation 
in terms of conventional matter models which are more 
familiar and intuitive compared with the "exotic" ones. 
A remarkable difference to scalar-field-driven inflation is 
the circumstance that our approach predicts an exponen- 
tial increase of the comoving entropy nsa 3 during the de 
Sitter phase. 

Our results may also shed new light on the old question 
whether or not a fluid bulk pressure may drive inflation. 
We recall that this issue has been discussed in the lit- 
erature from different points of view ( ||[35 38 , ) . 
While Pacher et al. || have shown that sufficiently high 
negative pressures cannot arise in a weakly interacting 
mixture of relativistic and nonrelativistic particles (see 
also jj7|), Lima et al. J35| pointed out that the situation 
may be different if the dilute-gas approximation is given 
up and causal thermodynamics is applied. Further in- 
vestigations along this line have confirmed the existence 
of inflationary solutions |36],[38]j|,||,|| , although there ex- 
ist general problems with their physical interpretation 
HHH, at least as long as cosmological particle produc- 
tion is not taken into account. The fact that a bulk pres- 
sure may phenomenologically represent certain quantum 
phenomena, especially particle production processes, is 
well known in the literature |39|-41|. As was remarked 
in |l3| , |l7|[ "conventional" bulk pressures, i.e. bulk pres- 
sures due to internal interactions, and effective bulk pres- 
sures resulting from particle production are separate ef- 
fects and both of them contribute to the overall dynamics 
of the system. 

We emphasize again that the quantity tt in the present 
paper is exclusively due an increase in the number of 
particles and not the "conventional" bulk viscous pres- 
sure of linear, irreversible thermodynamics, describing 
internal particle number preserving interactions. "Con- 
ventional" bulk pressures, generally equivalent to ,s =/= 0, 
have been excluded here by the assumption of collisional 
equilibrium. The possibility of a nonvanishing s due to 
the process of particle production was eliminated by the 
requirement of "generalized" equilibrium (see the discus- 
sion below Eq. (f2l|)). 

The above mentioned studies within the framework of 
causal thermodynamics relied on deviations from ther- 
modynamical equilibrium characterized by s ^ whereas 
the present considerations refer to a (generalized) equi- 
librium and imply s — 0. While there are limits for 
deviations from equilibrium in the mentioned nonequi- 
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librium appoaches (deviations up to second order) , there 
are no such restrictions in the context of this paper. In 
particular, there is no need of a requirement |7r| < p in 
our case to be well within the range of applicability of the 
theory as in conventional noncquilibrium thermodynam- 
ics. The present analogue of the bulk viscous pressure it 
which in conventional irreversible thermodynamics repre- 
sents (small) deviations from (collisional) equilibrium is a 
quantity without corresponding limitations. Instead it is 
determined by equilibrium conditions, equivalent to sym- 
metry requirements. This quantity 7r is directly related 
to the particle production rate which is traced back to a 
simple force on the (classical) microscopic level. It follows 
that 7r is completely determined by this force. The prob- 
lem whether a bulk pressure may drive inflation reduces 
to the question whether there exist microscopic forces 
on the particles, equivalent to a nongeodesic motion of 
the latter, which generate an appropriate macroscopic 
quantity tt. Our considerations show that a surprisingly 
simple force generates such a quantity. In this sense the 
question whether or not an effective bulk pressure may 
drive inflation is answered affirmatively. A final state- 
ment, however, requires the derivation of this force from 
an underlying quantum level which is beyond the scope 
of this paper. 

Having determined the selfinteracting force by the 
equilibrium conditions of the gas it is now also possi- 
ble to study the particle motion (|3^) explicitly. With the 
decomposition p % = Eu % + Xe 1 the left-hand side of the 
second equation ([||) may be written as 

Dp* _ dE_ ui + E Du l dX^ ^De l 
dr dr dr dr dr 

Contraction with u, yields 

Dp 1 dE , De 1 dE , ,D Ui 

Ui ~r~ = ~~r + Xu i~r~ = — Ae ~r~ ■ 

dr dr dr dr dr 

Taking into account that 

Du* , v n 



dr 



we obtain 



dE XE i . X 2 i n 



Dp 1 

~dV 



dr m 



Applying here the decomposition (|77j), the projected 
equation of motion 



Ui^f- = u t F l = -F 
dr 

may generally be written as 

A 2 „ A 2 



dE XE 
dr m 



3m 



-6 = F 



(78) 



(79) 



For homogeneous, isotropic universes with iii = <Xj n = 
the last equation reduces to 



dE_ 



3m 



-e = F . 



With dr = dt(m/E), X 2 = E 2 — rr, 
dE/dt = E, Eq. (|8(]) is equivalent to 



(80) 

= 3d/a and 



{E 2 - m 2 
E 2 — m- 



2m 



E 2 



(81) 



We discuss the last equation for three different cases: 
(i) F — 0, geodesic motion. We find 



A oc a 



(F = 0) 



(82) 



implying the expected behavior E cx a -1 for massless 
particles (photons) while the nonrelativistic energy e = 
E—m with e <C m of massive particles decays as e cx a~ 2 . 

(ii) F = (E — m) 0/3. The previously studied force 
with a linear dependence of F on the particle four- 
momenta is given by this expression for massive particles 
(m ^> T) in a homogeneous universe pif . In such a case, 
which implies the conformal Killing-vector property ( pd| ) 
and power-law inflation according to a oc t 4 / 3 |1, the 
solution of Eq. (pC|) is 



E 



m cx a 



F = (E - m) — 
K ' 3 



(83) 



The nonrelativistic energy of massive particles under gen- 
eralized equilibrium conditions in a (quasi-)linear force 
field decays linearly with the cosmic scale factor, i.e., in 
this case the selfinteracting force makes nonrelativistic 
particles behave like radiation. This may be regarded as 
the microscopic counterpart of the statement that radia- 
tion and nonrelativistic matter may be in equilibrium in 
the expanding universe, provided the number of matter 
particles increases at a specific rate . 

(iii) F = j-6. This is the case of interest here [cf. Eq. 



( p7| ) with Fi = and Fi from Eq. (62)]. It is obvious 
from Eq. j8C| ) that the force term on the right-hand side 
exactly compensates the second term on the left-hand 
side. Consequently, E vanishes, i.e. 



E = const 



3m 



(84) 



The self-interacting force prevents the particle energies 
from decaying with the expansion. Independently of the 
equations of state the particle energies are preserved in 
such a universe. 

With the result (|84|) we have completed the exact so- 
lution of our model of a quadratically self-interacting gas 
universe. It is the essential feature of this model that 
the same force which on the microscopic level makes the 
gas particles move at constant energy is responsable for 
an effective gravitational repulsion on the macroscopic 
level, implying an exponentially accelerated expansion of 
the universe. 
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V. CONCLUSIONS 

In this paper we introduced the concept of a 
"projector-conformal" timelike Killing vector (PCTKV) 
and discussed the corresponding fluid dynamics under 
the condition of minimal entropy production (generalized 
equilibrium). Such kind of equilibrium configuration re- 
quires a particle production rate which coincides with the 
fluid expansion rate. As a consequence the energy density 
of the fluid turned out to be stationary. A microscopic re- 
alization of this phcnomcnologically defined concept was 
given with the help of the kinetic theory for a classical 
gas in a force field. A quadratic (in the particle four- 
momenta) self-interaction of the microscopic gas parti- 
cles was shown to provide both the PCTKV property of 
Ui/T and "adiabatic" (or "isentropic" ) particle produc- 
tion. This force concept turned out to result in a com- 
prehensive picture of the gas dynamics both macroscop- 
ically and microscopically and allowed us to establish an 
exactly solvable model of a quadratically self-interacting 
gas universe. We found that generalized equilibrium un- 
der the conditions of spatial homogeneity and isotropy 
for such a configuration necessarily implies a de Sitter 
spacetime. We clarified in which sense an effective bulk 
pressure may drive exponential inflation. 
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